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We study the distribution of hydrogen and various hy- 
dride configurations in realistic models of a-Si:H for 
two different concentration generated via experimentally 
constrained molecular relaxation approach (ECMR) (TJ. 
The microstructure corresponding to low ( < 10%) and 
high (> 20%) concentration of H atoms are addressed 
and are compared to the experimental results with par- 
ticular emphasis on the size of H clusters and local envi- 
ronment of H atoms. 



The linewidths of the nuclear magnetic resonance 
(NMR) spectrum for the model configurations are cal- 
culated in order to compare to the experimental NMR 
data. Our study shows the presence of isolated hydrogen 
atoms, small and relatively large clusters with average 
proton-proton neighbor distance in the clusters around 
1.6 - 2.4 A that have been observed in multiple quantum 
NMR experiments. 
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1 Introduction Hydrogenated amorphous silicon (a- 
Si:H) is a very important material with a wide variety 
of technological applications ranging from photovoltaics 
to memory switching devices, thin film transistors, solar 
cells, optical scanners and numerous other electronic in- 
struments 0. The material has been studied extensively 
starting from its structural, electronic to vibrational and 
optical properties that produced a large volume of exper- 
imental data. In spite of this, our understanding of a-Si:H 
from theoretical point of view is very limited. Some of 
the very fundamental processes involving the material are 
poorly understood. For example, the creation of metastable 
defects in the material upon prolonged exposure to light 
by breaking of a Si-Si bond next to Si-H bond, the so- 
called Staebler-Wronski effect (SW) [4], is perhaps the 
most intriguing problem that is yet to be understood. Since 
many of the proposed microscopic mechanisms that at- 
tempt to explain the SW effect are based on bond-breaking 
models (6), it is crucial to understand the distribu- 
tion of H atoms and their dynamics in the amorphous 
network. Although there exists a number of experimental 
studies on hydrogen distribution in the material via nu- 



clear magnetic resonance (NMR) [7]|8]|9] and infrared 
(IR) spectroscopy (TO), there are very few theoretical stud- 
ies that address the problem explicitly flHl[T2l . 

In this paper we propose to calculate the hydrogen mi- 
crostructure of realistic models of a-Si:H at low and high 
concentration. Since one-dimensional NMR data cannot be 
uniquely mapped onto real- space distribution without the 
uncertainty of various fitting procedures, it is much more 
useful and convenient to employ first-principles models of 
a-Si:H to generate a distribution of hydrogen atoms. The 
NMR linewidths of the distribution can be obtained from a 
knowledge of the position of the spins in the network using 
a suitable approximation on the shape of the linespectra. 
This approach provides a direct route to study the number 
and the size of hydrogen clusters, their local environment, 
and the effect of the cluster size on the NMR linewidths. 
The dependence of H microstructure on the concentration 
of H atoms is also addressed and a comparison is made to 
experimental data. The plan of the paper is as follows. In 
section 2, we summarize the various experimental results 
from NMR experiments. Section 3 then describes briefly 
the a-Si:H models, and the resulting distribution of hy- 



Copyright line will be provided by the publisher 



2 



Timilsina et al.: Hydrogen microstructure of a-Si:H 



drogens at two different concentration with emphasize on 
sparse and clustered environment. In section 4, we calcu- 
late the NMR linewidths for the models and compared to 
experimental data, which is then followed by conclusion of 
our work. 



2 Summary of experimental results The two prin- 
cipal methods that are used to obtain experimental data on 
hydrogen distribution in a-Si:H are nuclear magnetic reso- 
nance (NMR) 00 and infrared spectroscopy (IR)COl. 
The former provides useful information about the local en- 
vironment of H atoms, whereas the latter plays a dominant 
role in determining the nature of various silicon-hydrogen 
bonding configurations. Experimental data suggest that 
H atoms can reside both in isolated and clustered phase 
and that the distribution can be highly inhomogeneous 
l9l fT3lL which depends on a number of factors such as 
preparation conditions, method of deposition, substrate 
temperature etc. A typical NMR spectrum of a device- 
quality sample shows the presence of both narrow and 
broad linewidths, which can be approximated as a con- 
volution of a Lorentzian (narrow) and Gaussian (broad) 
distribution. The dipolar interaction between the protons 
suggests that the narrow width is associated with a dilute 
or sparse distribution of spins, whereas a dense or clus- 
tered environment is needed to produce a broad linewidth 
in the resonance spectrum. Baum et al. demonstrated 
via multiple-quantum NMR experiments that for device- 
quality samples the clustering of H atoms increases with 
the increase of H concentration. The typical size of the 
clusters found in their work was of the order of 4-6 atoms, 
and with increasing concentration the clusters became 
physically closer to one another. NMR and IR measure- 
ments by Gleason et al. ifTH of a-Si:H prepared via plasma- 
enhanced chemical vapor deposition (PECVD) also con- 
firmed this observation. These authors reported the pres- 
ence of isolated, small and medium size clusters in the 
samples that gave rise to both broad and narrow linewidths 
of about 28-33 kHz and 6-12 kHz respectively in the con- 
centration range 10 - 16%. In other words, the picture that 
emerges from NMR and IR experiments on samples pre- 
pared by a variety of methods is that for device-quality 
samples H prefers to be distributed sparsely in the sample 
with occasional clustering at low temperature. While the 
clustering is more likely at high concentration, the mi- 
crostructure can be very inhomogeneous, and may consist 
of few large clusters even at very low concentration. An 
example of such inhomogeneous distribution of hydrogen 
atoms was reported by Wu et al. in hot-filament-assisted 
CVD deposited a-Si:H film [9]. They observed a signif- 
icantly large width of about 50 kHz in HW sample at 
very low concentration (2-3%) consisting of at least one 
H complex that was quite different from device-quality 
glow-discharged a-Si:H at 8-10%. This broad width is an 
indicative of the presence of large clusters in the sample 




Figure 1 (Color online) Hydrogen distribution and its lo- 
cal environment in a model of a-Si:H with 7% H atoms. 
The distribution of H atoms (in yellow/white) is sparse 
showing isolated, monohydride and dihydride configura- 
tions. The Si atoms are shown in blue/black with the 
exception of the one in the dihydride configuration (in 
green/white). 



that makes the distribution very inhomogeneous even at 
very low concentration. 

3 Microstructure from model a-Si:H In order to 
study hydrogen microstructure, we have used two realistic 
models of a-Si:H obtained via experimentally constrained 
molecular relaxation approach (ECMR)[ 1 ]. The method is 
based on a hybrid approach that has two important com- 
ponents: 1) experimental data can be directly utilized in 
building the model configurations, and 2) the model con- 
figurations are then driven via first-principles force-fields 
to obtain the minimum energy configuration. The method 
is unique in the sense that it is able to avoid some of the 
common pitfalls of both traditional first-principle molec- 
ular dynamics (small size, and large simulation time), 
and the uncertainty associated with constructing models 
solely based on limited experimental data (such as reverse 
Monte Carlo (2 |). We generate two models for our purpose 
of microstructure calculation. The first model consists of 
540 atoms with 7% H atoms, while the second one has 
611 atoms with 22% H atoms. The structural, electronic 
and vibrational properties of the models had been exten- 
sively studied by one of us and were compared to experi- 
mental data obtained from X-ray and Neutron diffraction 
measurements fT31 . Since the microstructure can be very 
sensitive to the details of the model construction, it is 
important to validate the properties of model networks 
with experimental data prior to addressing hydrogen mi- 
crostructure. 

The 540-atom model with 7% H atoms has very few 
defects and closely corresponds to a device-quality sam- 
ple at low concentration. The model has a clean electronic 
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band gap and the correct vibrational density of states show- 
ing the characteristic acoustic and optical peaks that have 
been discussed at length in Ref. fHll . A real-space analy- 
sis of this model shows the presence of both isolated and 
clustered phase of H atoms. In fig [I] we have shown the 
distribution of hydrogens in a region of linear dimension 
6 A that consists of 6-7 H atoms. The presence of isolated 
H atom, Si-H and SiH 2 configurations are also visible in 
the figure. An SiF^ configuration is explicitly shown with 
the Si atom in green/white. The local environment of this 
Si atom shows clustering of H atoms within the radius of 
3-4 A. The model consists of several such regions some 
of which have more dilute or dispersed atoms, while the 
others are relatively less sparse with small clusters of 4-5 
atoms. The proton-proton distance in this model is approx- 
imately lying between 1.6 A to 2.4 A. Our analysis shows 
that although there are few clustered phase (within 4-6 A), 
there are regions where the H atoms are separated by rel- 
atively large distances. The overall microstructure of the 
model is shown in fig. [3] without the embedding Si matrix 
(except silicon-hydrogen bonding) in order to highlight the 
sparse and densely occupied regions. It is evident that H 
distribution is quite sparse (except for few small clusters) 
consisting of few isolated H atoms, monohydride SiH and 
dihydride SiH 2 configurations. Gleason et al. ifTH studied 
a-Si:H films obtained by plasma enhanced CVD at vari- 
ous temperature in the concentration range 9 - 16 % of hy- 
drogen via multiple quantum NMR (MQ-NMR) and sug- 
gested the presence of small clusters in the samples. Baum 
et al. (8) also arrived at the similar conclusion studying 
device-quality samples via MQ-NMR in their works and 
reported the presence of cluster of 4 to 7 atoms. These ob- 
servations are consistent with the results obtained from our 
model network. The implication of disperse and clustered 
environment of H atoms in the network on NMR linewidths 
is discussed in the next section. 

We now consider the model with 22% hydrogen in the 
network. This value is somewhat high for the model to 
be considered as device-quality, nonetheless our intention 
is to study how the addition of more hydrogen changes 
the microstructure in the network. Experimentally it has 
been observed that with increasing concentration hydrogen 
atoms tend to form clusters until the concentration is high 
enough when small clusters are no longer visible lfT4l . In 
fig|2]we have shown the microstructure of a region of space 
having linear dimension of 6 A. The clustering of hydro- 
gen atoms surrounding an SiH4 configuration is explicitly 
shown in the figure. The local environment of this Si atom 
consists of 6-8 H atoms, and both monohydride and dihy- 
dride configurations are realized in this model. In addition 
to these, some isolated H atoms are also found to exist in 
the network. The entire region consists of about total 20 
H atoms showing both clustered and dispersed H atoms. 
Further analysis of different regions in the same model 
suggests that monohydrides can reside in both sparse and 
small clustered environment. The H-H neighbor distance 




Figure 2 (Color online) Hydrogen microstructure in a 
model of a-Si:H with 22% H atoms showing a cluster of 
H atoms (in yellow/white) consists of 6-10 atoms. A rep- 
resentative SiH 4 configuration is also shown with the Si 
atom in green/white. The other configurations present in 
the environment include isolated H atoms, SiH and SiH2. 



in this model is found to be in the range of 1.6 - 2.4 A and 
the distribution is shown in fig(5] A closer inspection of 
the local structure (with a cut-off 3.5 A) reveals that higher 
values at 2.1 A and 2.35 A mostly originate from SiH2 and 
SiH4 configurations. This can be easily understood con- 
sidering the bonding geometry of these configurations and 
taking into account of Si-H bond length and H-Si-H bond 
angle distributions. The microstructure of the entire model 
is presented in fig |4] showing hydrogen configurations that 
are realized in our model. This clearly demonstrates the in- 
homogeneous nature of the microstructure of the model. 

4 NMR linewidths of model a-Si:H In the preceding 
section we have discussed how the microstructure changes 
with increasing hydrogen content of the sample. Since the 
local environment of H atoms is reflected in the shape and 
the width of the NMR spectrum, one can obtain valuable 
information on microstructure from a knowledge of these 
quantities. However, direct calculation of NMR spectra is 
very difficult from first-principles, and one usually resorts 
to the approximation that the spectrum can be expressed 
as a convolution of Lorentzian and Gaussian distributions. 
NMR experiments on variously prepared samples at differ- 
ent concentrations reveal that at very dilute concentration 
randomly dispersed spins (via dipolar interaction) give rise 
to a narrow line width, whereas the presence of small clus- 
ters can be best attributed to the broader linewidth of the 
spectrum. These narrow and broad widths are generally as- 
sumed to be associated with the Lorentzian and Gaussian 
broadening of the spectrum respectively. The concentration 
of hydrogen in our model is sufficiently high to assume that 
the spectrum can be approximated by a Gaussian distribu- 
tion to estimate the broad linewidth. This can be confirmed 
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via calculation of r, the ratio of the fourth moment (^4) 
to the square of the second moment ({12) for the model net- 
works and compare to the same with a Gaussian distribu- 
tion. A value close to 3 (for ideal Gaussian) suggests that 
the spectrum can be approximated via Gaussian lineshape 
for which the FWHM is given by a = 2.36 021 • For 
an amorphous system with a random distribution of spins, 
the second and the fourth moments are given by lfT6l : 
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Tr[P x ] 
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where N is the number of spins, I = \ and Oju is the an- 
gle between the vector rjk with the applied magnetic field 
direction. A somewhat less accurate but simplified expres- 
sion for a is given by (9), 



a B (kHz) = 190(l£i 
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where is the distance between the spins in the unit of A. 
Similarly, an estimate of the narrow linewidth of the spec- 
trum coming from the dispersed H atoms can be obtained 
fromCSl, 
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where n is the number density of the spins without counting 
those in the clustered environment. 

In table I, we have listed the ratio of the moments, r 
and both the narrow and broad linewidths (FWHM) for the 
models for two different concentration. The average value 
of the moment is obtained from Eq|3]and is indicated in the 
table. The value of r is found to be close but less than 3 
in both the cases indicating the inhomogeneous nature of 
the distribution. This is expected in view of the inhomo- 
geneous microstructure that we have seen in figs[T]and[3] 
These values are consistent with those reported in the lit- 
erature (9). The broad linewidth associated with 22% 
H concentration is of the order of 50kHz, which is some- 
what high but is expected here because of the presence of a 
large cluster (cf. fig. [2]) of linear dimension 6.0 A. Our re- 
sults support the recent observation of a significantly large 
linewidth in hot-filament-assisted CVD deposited films by 
Wu et al. (9). On the basis of their NMR data, these au- 
thors suspected the presence of large H clusters in the sam- 
ple with few H nearest neighbors at a distance 1.6 - 1.8 
A. Together with the microstructure, one can infer that 



Figure 3 (Color online) The overall hydrogen microstruc- 
ture of the model with 7% H obtained from the ECMR 
method. The Si matrix is removed (except Si-H bonds) in 
order to display the nature of microstructure at low concen- 
tration. The Si and H atoms are shown in red (dark grey) 
and green (light grey) respectively. 



r 



Figure 4 (Color online) The overall hydrogen microstruc- 
ture of the model at high concentration (22% of H). Once 
again Si (in red/dark grey) and H (in green/light grey) 
atoms are displayed for clarity. Several H clusters are 
present along with few isolated H atoms, SiH and SiH 2 
configurations. Note that the microstructure consists of a 
combination of sparse and clustered environment that are 
observed experimentally. 



large H complexes (with 20 H atoms) are realizable, and 
when present they can produce a very broad linewidth in 
the spectrum. 

5 Conclusion The hydrogen microstructure in two 
realistic models of a-Si:H is studied at low and high con- 
centration of H atoms. The distribution of H atoms in both 
the models are found to be inhomogeneous that consists of 
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Figure 5 The distribution of H-H distances in a model of 
a-Si:H with 22% H atoms. A cut-off of 3.5 A is used to ob- 
tain the distribution. The peak at 1 .8 A is found to originate 
from the small clustered environment, whereas the peaks at 
2.1 and 2.35 A can be shown to appear from geometrical 
consideration of SiH2 and SiH4 configurations. 



few isolated H atoms, monohydrides and dihydrides that 
reside in both sparse and clustered environment. In addi- 
tion, few SiH4 configurations are found to realize in our 
first-principles model networks. At low concentration, the 
micro structure essentially consists of few small clusters of 
4-5 atoms and a sparse distribution of H configurations. For 
the network with high concentration of hydrogens, a clus- 
ter containing as many as 20 H atoms has been observed 
in a region of linear dimension 6 - 8 A. The theoretical 
linewidths obtained from the models using Gaussian and 
truncated Lorentzian broadening match with experimental 
data. In particular, our study confirms the presence of large 
clusters of hydrogen with 1.8 - 2.4 A as H-H neighbor dis- 
tances. The existing multiple quantum NMR data confirm 
our results. 
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Table 1 NMR linewidths for model a-Si:H 



Cone (%) 


Directions 


cjb (kHz) 


ajv (kHz) 


r 


7 


001 


16.34 (17.33) 


2.8 


2.8 




010 


19.68 (17.33) 


2.8 


2.6 




100 


17.80(17.33) 


2.8 


2.8 


22 


001 


52.38 (50.36) 


6.4 


2.9 




010 


52.15 (50.36) 


6.4 


2.9 




100 


49.92 (50.36) 


6.4 


2.9 
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